function f = parabola(y)
N = size(y,2);
f = zeros(N,1);
for i=1:N
    f(i) = 4*y(i)*(1-y(i));
end
end
